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1. TRANSFERENCIA DE CALOR POR CONVECCIÓN 


CAPA LÍMITE HIDRODINÁMICA O DE VELOCIDAD 


Sea una placa plana inmersa en una corriente fluida que discurre a velocidad u.. 
paralelamente a la placa. La viscosidad u del fluido da lugar a un esfuerzo cortante en 
las proximidades de la placa que provoca una reducción de la velocidad del flujo cerca 
de la misma. 


Se define como capa límite de velocidad o capa límite hidrodinámica a la región en la 
que existe gradiente de velocidad. El espesor 3 de la capa límite hidrodinámica es tal 
que el fluido experimenta una variación de velocidad que va desde O en la pared al 
99% de la velocidad de la corriente libre, definida como la zona que no está perturbada 
por la presencia de la placa: 


u u, Corriente libre 
u=0,99u, 
Capa límite 
hidrodinámica 
o de velocidad 


Borde de ataque 


Figura 1. Capa límite hidrodinámica. 


Con el incremento de la distancia al borde de ataque, el efecto de la fricción penetra 
tanto más en la corriente libre y el espesor de la capa límite hidrodinámica aumenta. El 
espacio fluido queda dividido en dos zonas: 


+ Corriente libre, llamada también como zona de flujo no viscoso o zona 
potencial: 
"= La velocidad de todas las partículas fluidas es igual a la de la 
corriente incidente. 
"=  Despreciable gradiente de velocidad (consecuencia de la anterior). 
"=  Despreciable esfuerzo cortante (consecuencia de la anterior). 


+ Capa límite o zona de flujo viscoso: 

"= Su espesor comienza siendo nulo en el borde de ataque. 

= Su espesor crece progresivamente con x. 

"= Su espesor en un punto es pequeño comparado con la coordenada 
del punto, esto es s(x)<<x. 

= La velocidad de las partículas fluidas es menor que la de la corriente 
libre y presenta cierta componente vertical. 

"= Elevado gradiente de velocidad. 

= Elevado esfuerzo cortante, que se hace máximo sobre la placa, por 
serlo el gradiente de velocidad. 
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Zonas de la capa límite 


Se distinguen dos zonas, una de flujo laminar y otra de flujo turbulento, cuyas 
características más significativas son: 


Vo Vo [o.0) 
3 1 
LS Sl Núcleo 
LS Mea turbulento Cc 
EA . apa 
2] a, 3 límite 
Es E 

A IA Subcapa 
GD + laminar 


Laminar Turbulento 
Figura 2. Zonas de la capa límite 


Laminar: 


"= Se extiende desde el borde de ataque hasta el denominado punto de 


transición, lo que matemáticamente se expresa como 
YX ek 
Re, = A Re. 
17, v 
=  Lavelocidad de las partículas coincide con su velocidad media, v = v. 
4,91 x 


= Espesor de la cala límite es: Y = 72 


Xx 


Re 


Turbulento: 


"=  Seextiende desde el final de la capa límite laminar (punto de transición), hasta 
el borde de salida de la placa. Matemáticamente Re, > Re. 

=  Lavelocidad tiene una componente aleatoria v' o de vibración espacial que se 
superpone a la velocidad media V, esto es, v=V+V". 

=  Esel tipo de flujo que existe en la mayoría de las aplicaciones de interés. 

= Se distingue una zona próxima a la pared con características muy parecidas a 
las de la capa límite laminar y a la que se denomina subcapa laminar. Dentro 
de ésta las partículas carecen de vibración (v'=0), en contraste al núcleo 
turbulento, donde sí existe tal vibración. 

0,377 x 

Re'* 
29x 


Re 


x 


= Espesor de la capa límite: d = 


= Espesor de la subcapa laminar: J, = 


La transición de flujo laminar a turbulento tiene lugar en un punto de la placa tal que 
Re, alcanza un valor crítico Re¿, comprendido en el intervalo 5-10” a 6-10*%, según la 
mayoría de los autores. El que Re¿ tome un valor más próximo a uno u otro extremo va 
a depender del nivel de perturbación de la corriente, de manera que cuanto mayor sea 
dicha perturbación, Re, estará más cerca de 5-10", y viceversa, Re, se aproximará a 
6-10* cuando el nivel de perturbación sea tanto menor. Esta perturbación puede ser 
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inducida bien por la placa, mediante vibración o introduciendo rugosidad superficial, o 
bien por interposición de pequeños elementos que generan vorticidad, como un fino 
hilo dispuesto perpendicularmente a la corriente incidente. 


La determinación de la coordenada x. en la que tiene lugar la transición de flujo 
laminar en turbulento resulta de: 


xv Re v 
Re, => x == 
v v 


00 


De donde los valores extremos mínimo y máximo de x.: 


Re. v 10%y Re. v 10 
= X = 


max Cmin 
v v v v 


00 00 00 00 


Que expresado de una forma gráfica: 


Placa de baja rugosidad 
Placa no vibra 
Corriente incidente no perturbada 


C Cmax 
Vórtices 
corriente ; 
incidente Rugosidad 
superficial e 
6 S Vibración de 
C 


Placa de alta rugosidad 
Vibración placa 
Corriente incidente perturbada 


XIX 


4 min a 


Laminar Turbulento 


Figura 3. Transición de la capa límite laminar a turbulenta. 


CAPA LÍMITE TÉRMICA 


Si la temperatura T., de la corriente incidente es distinta de la temperatura T, de la 
placa, en las proximidades de ésta se desarrollará una zona de transición entre ambas 
temperaturas. 


Se define como capa límite térmica a la región en la que existe gradiente de 
temperatura. El espesor 6, de la capa límite térmica es tal que el fluido experimenta 
una diferencia de temperatura igual al 99% de la diferencia existente entre la placa y la 
corriente libre, definida como la zona que no está perturbada térmicamente por la 
presencia de la placa: 


Corriente libre 
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= 0,99 
Capa límite 
A: térmica 
AN Ly) 
ELILA AA AA AA AAA ANAL ANA LALALA LAA, 
q” q” 
TL LS, 
L 


Figura 4. Capa límite térmica 


Con el incremento de la distancia al borde de ataque, el efecto de la transferencia de 
calor penetra tanto más en la corriente libre y el espesor de la capa límite térmica 
crece. 


Características de la corriente libre: 

=  Despreciable gradiente térmico 

=  Despreciable transmisión de calor 
Características de la capa límite: 

= Elevado gradiente térmico 


= Elevada transferencia de calor 


Indicar que la existencia de zona laminar y turbulenta sobre la capa límite 
hidrodinámica, también da lugar a la misma división sobre la capa límite térmica. 


Capa límite 
térmica 


Laminar Turbulento 
| taminar | uranio 


Respecto a los espesores de las capas límite, decir que el de la capa límite 
hidrodinámica está directamente relacionado con la difusividad cinemática O 
viscosidad cinemática v=u/p (u: viscosidad dinámica, p: densidad), de tal manera que 
cuanto mayor es v, tanto más elevado es dicho espesor. De esta forma, la capa límite 
de la figura a es representativa de una v más pequeña que la de la figura b, supuesto 
que en ambos casos la velocidad de la corriente libre es la misma. 
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a 


Figura a Figura b 


Idéntica relación existe entre la difusividad térmica «=k/pc, (k: conductividad, p: 


densidad, c,: calor específico) y el espesor de la capa límite térmica, es decir, el 
espesor de la capa límite térmica crece con a. En cualquier caso, las capas límite 
hidrodinámica y térmica son delgadas comparadas con las dimensiones de cuerpo 
sobre las que se desarrollan. 


El número de Prandtl (P,=v/a) es una medida relativa de los espesores de las capas 
límite hidrodinámica y térmica. En concreto, para una placa plana calentada a 
temperatura uniforme T;,, la relación entre los espesores de las capas límite 
hidrodinámica d y térmica 8;: 


Lo 3/p 
Ó, j 
pudiéndose generalizar que: 
si P>>1> v>>a > 0>>, 
si P¡<<L> v<<a => 0<<6; 
y, por tanto: 
Metales líquidos (P,-0,001): 9/85;-0,1 
Aire (P,-0,7): 8/81 
Agua (P,-5): 8/8;-1,8 
Aceites lubricación (P,-1000): S/8;-10 


Variación cualitativa de h(x). El flujo de calor (potencia por unidad de área) en un 
punto cualquiera de la placa se puede expresar como: 


q =h(T, —T,) (1.1) 


donde h=h(x) puede variar a lo largo de la placa. Dado que el fluido en contacto con la 
placa está en reposo, la transferencia de calor a través de esta capa fluida tiene lugar 
exclusivamente por conducción, pudiéndose expresar mediante la ley de Fourier: 
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(1.2) 


y=0 
siendo k; la conductividad del fluido. Combinando las ecuaciones (1) y (2), h queda 
como: 


k 
Ms. 


— 1.3 
T,_-—T, 0y si 


y=0 


Una aproximación del gradiente de temperatura en la pared se puede obtener 
sustituyendo la ley de temperatura en la capa límite por una relación lineal (línea roja 
en la figura 4): 


rl _E=-T 
(8) 
Wly.o 0% 
y por tanto: 
k az k 
e e: al 


Te En Lo Ó, t 


Y dado que 3; aumenta con la distancia x, h(x) disminuye con x, lo que gráficamente se 


recoge en la siguiente figura: 
/ hCg « qx) 


Laminar Turbulento 


Figura 6. Evolución de h y 8; a lo largo de una placa plana. 


donde puede apreciarse que el coeficiente de convección experimenta un salto al 
pasar de flujo laminar a turbulento. Esto justifica el hecho que en el flujo turbulento, ya 
sea sobre una placa plana o sobre cualquier otro cuerpo, la potencia transferida por 
unidad de área o flujo de calor (q”) es mayor que en flujo laminar. 


Por otro lado, combinando (2) y (3”): 


que indica que q”(x) disminuye a lo largo de la placa de una forma cualitativa 
semejante a la de h(x), experimentando también un salto discreto al pasar de la zona 
laminar a la turbulenta. 
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Definición del valor medio de h. La potencia calorífica total evacuada por un tramo 
de placa rectangular de longitud L y anchura b: 


q=5,q'dA=| h(T, —T,,)dA 
Para el caso en que la temperatura superficial T¿ sea uniforme en la placa: 
q= |, ACT, —T.) dA= (T, =cte)=(T, -T,)| hdA=hA(T, —T,) (1.4) 
de donde: 
h = e fr dA T.=cte 


O bien, teniendo en cuenta que A=b-L, siendo b el ancho de la placa y L la longitud de 
la misma: 


h=-= y hdx T.=cte (1.5) 


A pesar de que esta expresión de h han sido obtenida bajo el supuesto que la 
temperatura superficial es constante (T.=cte), se considera una expresión general, 
válida tanto si dicha temperatura superficial es constante como si no lo es. 


Potencia calorífica total intercambiada. La definición de h hace posible el cálculo de 
la potencia térmica total que intercambia la placa plana de una forma sencilla 
mediante: 


q= hA(T, -T,) T.=cte (1.6a) 


la cual solamente es válida para T.=cte. Cuando esta condición no se cumple, la 
expresión de q toma la forma: 


q=hA(T, -T,) T.x*cte (1.6b) 
donde T; —T,, representa el valor medio de la diferencia de la temperatura entre la 
superficie y el fluido. 


Definición del Número de Nusselt. Partiendo de la ecuación (1.3): 


k 
h(x) = a 


0 1.3 
dd AS a 10 es 


0 


Se va adimensionalizar el término or /0y|,.. para lo cual se toman como variables 


características la longitud L de la placa y T.-T..: 


y =2 >  dy=Ldy 
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T" = E ES dT =(T. -T )dT' 
TE == ( S | 
y, por tanto: 
k * 
o (1.7) 
L 0y Les 


donde la magnitud adimensional - or” /ay| . ¿Se denomina número de Nusselt (Nu), 
7 


y tiene el significado del gradiente adimensional de temperatura del fluido en contacto 
con la pared. Haciendo uso de la definición de Nu, la ecuación (7) se puede reordenar 
como: 


L 
Nu = NN X hL_Kk da y RoNDUCTIVAFLUIDO (1.9) 
Oy E k; E ReoNvecrivA 
hA 


de donde se obtiene la interpretación física de Nu, entendido como el cociente entre la 
resistencia térmica a la conducción del fluido y la resistencia convectiva. 


Al igual que se ha hablado de coeficientes de convección local h(x) y medio h, se 
puede hacer lo propio con el número de Nusselt. Nusselt local (Nu(x) o Nu,) es el que 
toma en un punto determinado y Nusselt medio (Nu) el valor promedio en toda la 
superficie considerada: 


Ñ Ah) x 


Nusselt local: |Nu(x) = Nu, A 


(1.10) 


f 


Nusselt medio: | Nu = — (1.11) 


La relación entre ambos se obtiene partiendo de la ecuación (5), sustituyendo h y 
h por los valores dados por (10) y (11): 


Nu k; A [ Nu, k; = L Nu, 


d. Nu = dx 1.12 
E xo > Í (1.12) 


0 Xx 


La ecuación (12) tiene la utilidad del cálculo de Nua partir del conocimiento de Nu,. 


Proceso de cálculo de la potencia térmica en problemas de transferencia de 
calor por convección. Las etapas seguidas suelen ser las siguientes: 


1. Determinación del número de Nusselt medio Nuasociado al problema en 
cuestión (vendrá condicionado por la geometría y la hidrodinámica). 
2. Cálculo del coeficiente de convección medio (ecuación 11 reordenada): 
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k; N, 
L 
3. Determinación de la potencia térmica propiamente dicha (ecuaciones 6): 


h= 


q=hA(T,-T,) 6  q=hA(T,-T,) 


Dependencias funcionales de las correlaciones en convección. Lo dicho aquí está 
relacionado con el apartado a del punto anterior. Según que la convección sea 
forzada, natural o mixta las dependencias funcionales de Nu cambian: 


Convección Forzada: Nu= f (geometría, Re, Pr) 
Convección natural: Nu= f (geometria, Gr, Pr) 


Convección Mixta: Nu = f (geometría, Re, Pr, Gr) 


donde Re, Pr y Gr son los denominados números adimensionales de Reynolds, 
Prandtl y Grashof, cuyos significados físicos son: 


u,L A 

Re = Fuerza de inercia entre fuerza viscosa 
12) 
T, -T,JL 
Gr = 9PE, TIE Fuerza de empuje entre fuerza viscosa 
12) 

pe Relación entre las difusividades cinemática y 

Q térmica 


Desde un punto de vista conceptual las correlaciones en convección forzada y natural 
pueden expresarse como Nu=CRe” Pr” y Nu=C'Gr” Pr” donde C y C' son 
cantidades constantes que dependen de la geometría del cuerpo (superficie plana, 
flujo interior a un conducto, flujo exterior a un conducto, etc). Los exponentes m, n, m' 
y n' pueden determinarse matemáticamente mediante el estudio exhaustivo de la capa 
límite, pero esto escapa del objeto de este curso, aunque en cualquier caso pueden 
considerarse conocidos. Por tanto, la determinación de la única constante (C Ó C”) 
requeriría de un solo ensayo. A pesar de lo dicho, estas ecuaciones sólo deben 
utilizarse para analizar la tendencia de Nusselt con los números adimensionales 


implicados y no para el cálculo preciso de Nu de cara a la resolución de problemas de 
transferencia de calor. En su lugar deben emplearse algunas de las múltiples 
ecuaciones de ajuste a datos experimentales, denominadas correlaciones, cuya 
tendencia es la misma que la dada por Nu=CRe” Pr" o Nu=C'Gr” Pr”, y que 
será presentadas al final de cada capítulo de convección. 
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CORRELACIONES PLACA PLANA EN CONVECCIÓN FORZADA 


Re, == Re, === Re, == Re¿=5:10*-6-10* 


T_ +T 


00 Ss 


Propiedades a temperatura de película: T, = 


TEMPERATURA SUPERFICIAL UNIFORME (T.=cte) 
ES e AE 
Potencia calorífica: q =h -A-(T, —T,) h = 7 Jr dx 
0 


Flujo laminar (0O<x<L si Re,<Re.): 


Nu, = 0,332 Re!” Pr'* Pr>0,6 
Nu =0,664Re'” Pr'* Pr>0,6 


Churchil €. Ozoe: 


0,3387 Re”? Pr 
Nu, = z Re,Pr>100  VPr 


0,0468 Y" 
1+ C ] 
Pr 
Mu 0,6774 Re”? pr 
0,0468 4" 
1 > ] 
Pr 
Metales líquidos: 


Nu, =0,565(Re, Pr)" Re,Pr>100  Pr<0,05 
Nu =1130(Re, Pr)” Re, Pr>100  Pr<0,05 


Re_Pr>100  VPr 
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Placa con extremo no calentado (x>é): 


Us Tos 1/2 p, 13 
0,332 Re” Pr 
zz Nu, =" z Pr>0,6 
a ga 43 
5 ó 
¡E AA 
(7 
Flujo turbulento (x.<x<L): 
Nu, =0,0296Re*" Pr'* 0,6<Pr<60 
Nu=0,037(Re?”—Re**)Pr' 0,6<Pr<60 
Placa con extremo no calentado (x>g): 
_ 0,0296 Re** Pr!" 
a 9/10 112 
1 3 
X 
Flujo mixto (O<x<L con tal que Re¿<Re;,): 
Chilton-Colburn: Nu = [0,664Re"+0,037(Re?*-— Re” J]Pr”* 0,6<Pr<60 
Re,<10* 
R 4/5 
Witaker: Nu=0,664Re”” Pr*+0,036Re? Pr** E - e) 0,7<Pr<400 
e, 
Re,<3-10” 
FLUJO DE CALOR SUPERFICIAL UNIFORME (q”=cte) 
: e . O a le 
Potencia calorífica: q=q A=h Alr, E) h = ¿ro dx 
0 
Flujo laminar (O<x<L si Re,<Re,): 
Nu, =0,453Re'” Pr” Pr>0,6 
Nu =0,906 Re!” Pr'* Pr>0,6 
Flujo turbulento (x.<x<L): 
Nu, =0,0308Re?*” Pr'” 0,6<Pr<60 
Nu=0,0385(Re?*- Re?” )pr”* 0,6<Pr<60 


Flujo mixto (0<x<L con tal que Re¿<Re;): 


Nu =[0,906 Re”? +0,0385(Re?*- Re** )]pr'*  0,6<Pr<60 
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FLUJO TANSVERSAL FORZADO A UN CILINDRO ISOTERMO 


E a TA 
Propiedades a temperatura de película: — ] z 


Churchill €. Bernstein (Pr>0,5): 


0,62 Re”? Pr? 
2/3 y4 
(35) 
Pr 
E 0,62 Re”? Pr'P de 
Nu =0,3+ 1 as ] 


041% ys 282000 
14+| = 
E) 


0,62 Re"? pr Re y 


+ 
041% ye 282000 
1+| 2 
e 


Nu= z 
0,8237 — In(Re Pr)” 


Nu=0,3+ 
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Nu=0,3+ 


Nakai €. Okazaki: 


Re-Pr<0, 2 


Re<10* 


2-10%<Re<4-10* 


4-10"<Re<5-100 


Hilpert: Nu=CRe” Pr” 0,4<Re<4-10* Pr>0,6 
Re C m Re Cc m 

0,4-4 0,989 | 0,330 Ye, <> fD | 5.000-100.000 | 0,246 | 0,588 
4-40 0,911 | 0,385 Dto | 5.000-100.000 | 0,102 | 0,675 
5.000-19.500 | 0,160 | 0,638 
40-4.000 0,683 | 0,466 O 1D | 19500-100.000 (0,0385 | 0.782 
4.000-40.000 [0,193 | 0,618 O to | 5.000-100.000 | 0,153 | 0,638 
40.000-400.000 | 0,027 | 0,805 D fo | 4.000-15.000 | 0,228 | 0,731 


Cilindro circular con flujo cruzado. 


Cilindro no circular con flujo cruzado. 
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Zhukauskas. Cilindro circular. Propiedades a T., excepto Pr, a Ts. 


ya 
Nu=CRe” ES 0,7<Pr<500  1<Re<10* 
r 


S 


Pr<10 > n=0,37 Re [e m 
Pr>10 > n=0,36 1-40 0,75 0,4 
40-1.000 0,51 0,5 


1.000-100.000 0,26 0,6 
200.000-1.000.000 | 0,076 | 0,7 


FLUJO TANSVERSAL FORZADO A UNA ESFERA ISOTERMA 


D 
Ts Re = Ma 
v 
Veo 
To == hD 
No TS 
k 
Propiedades a T., excepto us a Ts. 
1ja 3,5< Re < 7,610* 
Nu =2+(0,4Re'”+0,06 Repo Le 0,71 < Pr < 380 
Us 
¡A 
Us 
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2. CONVECCIÓN NATURAL O LIBRE 


Es debida a la acción conjunta de la diferencia de densidad en el fluido y de la 
gravedad. Las diferencias de densidad entre puntos de una masa fluida pueden ser 
debidas a gradientes de temperatura o bien a gradientes en la composición. Esto 
último es lo que justifica, por ejemplo, la elevación del aire húmedo debido a la menor 
densidad del vapor de agua. El gradiente de densidad de origen térmico en un fluido 
se cuantifica a través del coeficiente de dilatación volumétrica: 


__1/(0p 
pa E), 


$ es una propiedad de las sustancias que en el caso de los gases ideales, y cuando no 
se encuentre tabulado, se puede determinar con buena aproximación mediante: 


p= | 2 
pxoT),  poOrLRT), pil RT?) T 


siendo T la temperatura absoluta; para los líquidos f se encuentra tabulado. 


Los flujos de convección natural pueden ser externos o internos, dependiendo de si 
están confinados parcial o totalmente por superficies. Las velocidades habituales 
asociadas a la convección natural son relativamente pequeñas, máximo en torno a 
unos 2 m/s, por lo que el coeficiente de convección natural suele ser más pequeño 
(-10 W/m*K en gases y -50 W/m*K en líquidos) que en convección forzada. Como 
consecuencia de esto el espesor de la capa límite térmica, en términos generales, es 
mayor en convección natural que en convección forzada para las velocidades 


habituales de la corriente. 
Fluido 
estático 
Too» Poo 
é 
AAN y 


Convección natural Convección forzada 


Figura 1. Capas límite sobre una placa plana en convección natural y forzada. 


Se demuestra que la velocidad característica que se desarrolla en convección natural 
sobre una superficie plana como la representada en la figura 1 viene dada por 


u=./B AT gL , lo que pone de manifiesto que aumenta con el coeficiente de dilatación 


volumétrica, con la diferencia de temperatura entre el fluido y la superficie y con la 
gravedad (en ausencia de gravedad no hay convección natural). 
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Con relación la figura 1, que representa las capas límite hidrodinámicas en sendas 
parades verticales, pueden hacerse las siguientes consideraciones: 


= Mientras que en convección forzada la velocidad del fluido es nula en la pared 
y u,, fuera de la capa límite, en convección natural la velocidad del flujo es nula 
en la pared y en los límites de la capa límite hidrodinámica. 


= El espesor de la capa límite hidrodinámica es tanto mayor cuanto menor es la 
velocidad característica del fluido (inducida naturalmente en convección natural 
o de manera externa en convección forzada). Dado que dichas velocidades 
suelen ser menores en convección natural, el espesor de la capa límite suele 
ser más grande en ésta última. Esto mismo rige para la capa límite térmica, y 
dado que el flujo de potencia térmica (q”) a través de la pared aumenta con 
espesores de capa límite térmica tanto más pequeños (ver capítulo 
introductorio de la convección), se deduce que en convección forzada la 
potencia térmica por convección suele ser mayor que en convección natural, o 
lo que es lo mismo, el coeficiente de convección es mayor en convección 
forzada que en convección natural. 


Al igual que en el flujo forzada paralelo a una placa plana podía tener lugar una zona 
laminar y otra turbulenta, en flujo natural también puede darse dicha circunstancia. 


Flujo turbulento 
h¿ = Constante 


Transición - 


Da, lo 


4 


Distancia desde la base de la pared, x 


Lamina —————=]-——— Turbulento 


Flujo laminar 
he ox =1M 


+ 
. 
PARITARIA TUTTI VTA 


po 


Pared caliente 


Figura 2. Capas límites laminar y turbulenta sobre puna placa plana en convección natural. 


No obstante, mientras en el flujo turbulento forzado h(x) disminuía al alejarse del borde 
de ataque, en el flujo natural turbulento h(x) permanece prácticamente constante a lo 
largo de la placa. 


Cuando se dan las condiciones apropiadas, el calentamiento de un fluido suele 
producir una corriente natural ascendente, o una corriente descendente si el fluido se 
enfría. Esto es sí dado que la mayoría de los fluidos experimentan una reducción de la 
densidad al calentarlos a presión constante y viceversa, esto es, un aumento de la 
densidad cuando se enfrían. Hay, no obstante algunas excepciones a esta regla. La 
más característica es el agua entre OC y 4*C, intervalo en el presenta un 
comportamiento inverso al descrito, pues en este rango al aumentar la temperatura la 
densidad aumenta en vez de disminuir, lo que motiva un movimiento natural inverso al 
descrito durante los procesos de enfriamiento y calentamiento. 


El número de Reynolds no aparece explícitamente en la convección natural puesto 
que la velocidad inducida o velocidad característica del flujo no suele ser conocida a 
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B-AT-g-L 
v? 
Rayleigh Ra = Gr -Pr, siendo $ el coeficiente de dilatación volumétrica, AT la diferencia 


de temperatura característica, y la gravedad, a la difusividad térmica y Pr el número de 
Prandtl. 


priori. En su lugar interviene el número de Grashof, Gr = o el número de 


Flujo natural externo 


La figura 2 muestra el desarrollo de la capa límite en una pared vertical a distinta 
temperatura que la corriente libre y la evolución del coeficiente de convección h, a lo 
largo de una pared calentada. En tal configuración, la transición de la capa límite 
laminar a turbulenta ocurre cuando Ra,-10*. 


El flujo de calor en convección natural externa tiene la misma forma que en flujo 
forzado externo q =h(T, —T,,), con la única diferencia que en convección natural la 


velocidad no es inducida externamente sino por el propio fenómeno de tranferencia de 
calor. En la expresión anterior Ts y T., son las temperaturas de la pared y de la 
corriente libre exterior. 


En flujo externo las correlaciones entre Nu y Ra tienen la forma cualitativa: 


Nu=C (Gr - Pr)” 
O bien: 
Nu=C Ra” 
donde n=1/4 y 1/3 para flujos laminar y turbulento, respectivamente y donde Ra es el 


número de Rayleigh. Es una consecuencia inmediata de la ecuación anterior que en 
flujo turbulento, h es independiente de la longitud característica sobre una placa plana: 


h 


1 1 
Nuk  C(GrPr*k -AT-q-13_ Y AT - 3 
E y _C(GrPrr"k_C B:AT-g lp, k=C B:AT-9 p, el 
L L L Vv V 


Flujo natural interno 


Es el que tiene lugar en un recinto cerrado en el que existe diferencia de temperatura 
entre las superficies que lo limitan (paredes con cámara de aire, colector solar de 
placa plana, ventanas con doble acristalamiento, etc). El flujo de calor se expresa a 
partir de las temperaturas características de las paredes caliente (T;, hot) y fría (Tc, 
cold): 


po id 


El parámetro de longitud característica para el cálculo del número de Nusselt es la 
separación L entre las placas. Si la diferencia de temperaturas Ty-T es menor que el 
valor crítico requerido para que el flujo se vuelva inestable, el calor se transfiere por 


conducción y, por tanto, q =, —T,), de donde h=k/L, y el número de Nusselt 
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Nu= TE =1, valor que constituye el límite inferior del número de Nusselt en las 


correlaciones de convección natural entre dos paredes planas. 


La inestabilidad del flujo interno a partir del cual se desarrolla la convección natural 
está dado por el número de Rayleigh y depende de la geometría del recinto, siendo 


característicos los siguientes valores: 


Paredes planas horizontales de espesor L de hueco: 


AS === 3 
Ra 700 === 


Paredes planas verticales de espesor de hueco L: 


Ra,=1000 


Espacio entre cilindros y esferas concéntricos: Ra, , =100 
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CORRELACIONES CONVECCIÓN NATURAL 
Superficies a temperatura T, uniforme 


Cilindros y esferas de diámetro D: 


_BATGD", BATgGD ELE: 


Ra S 
v va D 


Paredes y placas de longitud L y recintos internos con separación L entre caras: 


3 3 WN, 
a a PAI pen y Nu-k 
v v-a L 
FLUJO EXTERNO 
T.¿T. 
Propiedades evaluadas a la temperatura de película: T, == ] 
Pared vertical 
Longitud característica L: altura de la pared 
Fluido 
estático 
pl 
Churchil € Chu: 
= aj 
Nu =0,68+ pe 5 Ra<10? Laminar 
Ln” 
14| 222 
Pr 
2 
ER ys 
Nu =+0,825+ 1 ab7 Ha Ra>10? Turbulento 


8/27 
| | 
1+ 
Pr 
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Pared inclinada 
Longitud característica L: longitud de la pared en el sentido del flujo. 


Considerar: Ra, = Ra-cosO' = Ra «sin 0 


Figura 4,35 Patrones de flujo alrededor de placas inclinadas (a) calentadas y (b) enfriadas. 


0,670 Ra!* 


Nu =0,68+ Ras<10* Laminar 
0,499 1% 
ey" Pe 
2 
E TRA 
Nu =30,825 + DAOT BOS Ras>10* Turbulento 
04920" 
1 At ES | Omin=30% 
Pr 


Estas correlaciones solamente son válidas para las caras A indicadas. Para las caras 
B hay que consultar bibliografía específica. 


Flujo alrededor de un cilindro horizontal de gran longitud 


Longitud característica D: diámetro del cilindro 


Morgan: Nu =C-Ra" 
Ra Cc n 
10% a 10? 0,675 0,058 
10% a 10* 1,02 0,148 
10% a 10* 0,85 0,188 
10% a 10” 0,48 0,25 
10” a 10* 0,125 0,333 
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0,387Ra"* 


Churchill €, Chu: Nu, =30,60+ == Ra<10** 
| 0559 y" 
1+| 
Cr) 

Flujo alrededor de una esfera 

Longitud característica D: diámetro de la esfera 

a 0,589Ra** al 
Churchill: Nu=2+ Ra<10 Pr>0,7 


9/16 149 
| e ES | 
Pr 


Superficies horizontales 


Y ai. 
e 
Fluido, 74 


Placa, 7, 
> Dl 
Fluido, 1) 
(a) (b) 
Fluido, 7.., 
Placa, T, / 


CARAS 


WAY 


Fluido, 7... 


a) Fluido por encima de superficie con fuerza de flotación descendente, b) Fluido por 
debajo de superficie con fuerza de flotación descendente, c) Fluido por encima de 
superficie con fuerza de flotación ascendente, d) Fuido por debajo de superficie con 
fuerza de flotación ascendente. 


Longitud característica L=A/p, donde A es al área de la superficie y p su perímetro. 
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Fluido por encima de superficie con fuerza de flotación descendente y Fuido por 
debajo de superficie con fuerza de flotación ascendente 


Nu =0,27Ra”* 10*<Ra<10% 


Fluido por debajo de superficie con fuerza de flotación descendente y Fluido por 
encima de superficie con fuerza de flotación ascendente 


Nu =0,54Ra”* 10'<Ra<10” 
Nu=0,15Ra”? 10'<Ra<10* 


Objetos de forma arbitraria 


Lienhard: Nu = 0,52R,** donde debe emplearse como longitud característica la de 
la | capa límite térmica (L=1D/2 para cilindros y esferas). 
FLUJO INTERNO 


T, +T, E E 
Propiedades evaluadas a la temperatura: a (Tu: superficie caliente, Tc: 
superficie fría). 


La longitud característica L viene dada por el espesor del hueco entre paredes. 


Recinto hueco rectangular inclinado (H/L>10) 


lado caliente 


0<0<60* ; O<Ra<10? 


. 1,6 13 
Nu mtozaa1 Pf, Malta) |, (Ens) 1 


a-cosO Ra(cos0) 5830 


donde los términos entre paréntesis deben hacerse nulos si resultan negativos. 


0 = 60"; O<Ra<10” 
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Nu, = max(Nu,, Nu, | 


0,0936Ra””* 


0,5 
1+ 5 
l a (Ra/3160)** |” | 


Nu, =| 0,104+ 04175 | pgozs 
H/L 


60%< 9 < 90” 


0 =90"; 10é<Ra<10” 
Nu,, =max(Nu,, Nu,, Nu,) 


Nu, =0,0605Ra” 


3113 


0,104Ra””* 


R 0,272 
Nu, =0,242| Y 
H/L 


Nu¿y =1 para Ra<1000 


Recinto hueco vertical (superficies horizontales aisladas) 


aislado 
7 


Pr 0,29 
Nu =0,8 Ra 
0,2 + Pr 
0,28 1/4 
Nu=0l "RE 
,2+Pr L 


0,3 
Nu =0,42Ra”* pro? (E 


Cilindros largos concéntricos 


Cilindro interior 
D;,T; 


Cilindro exterior 


< HE <2 


10? <Pr<10? 
107 < Ra - Pr 


0,2 + Pr 


2<H/L<10 
Pr<10? 
10% < Ra<10*” 


10<H/L<40 
1<Pr<2-10* 
10* <Ra<10” 


<100 


Do, To 
Patrón de flujo 
PT 
: BEN D,-D. 
Longitud característica: L =—% E 
q A 
2D 
In 
D, 
27 k,¿L 
k,¿ =k Ra” 
k le Veas Ñ 
« =0,368 ———| (Raj 10?%<Ra'<10” 
k 0,861 + Pr 
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Esferas concéntricas 


; a D,-D, 
Longitud característica: L == 
AT 
E 
D, D, k¿ 7 


ya 
A ] (Ra")” 


0,861 + Pr 
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Patrón de flujo 
SAL 


10? < Ra" <10* 


3. CONVECCIÓN MIXTA 


Es la que tiene lugar cuando la convección natural y la forzada están presentes 
deforma significativa ambas. Algunos esquemas de convección mixta, entre los 
muchos posibles, se muestran en la tabla siguiente: 


T. 
ñ “Te O 


poa paz 


(a) Placa vertical (b) Cilindro (c) Esfera 


Convección natural y forzada a favor: T¿>T, y v? Ó T,¿<T, y wY 


(d) Placa vertical (e) Cilindro (/) Esfera 
Convección natural y forzada en oposición: T¿>T. y vi Ó T.<T. y v? 


vT, VT, 
es ————> 
E O O A 
£] £ | 


(g) Convección natural y forzada compuestas transversalmente: 
ToTe y VO Ó Te<Te y VO 


(¿) Cilindro * (7) Esfera (k) Cilindro 


Convección natural y forzada compuestas transversalmente: 
Tse>Te y VO Ó Ts<Te y VO 
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Si Nu, representa el número de Nusselt de la convección natural y Nu el Nusselt de 
la convección forzada, Churchill propone las siguientes correlaciones para la 
determinación del Nusselt mixto ( Nu ): 


1. Flujo de capas límite laminares o turbulentas con convección natural y forzada 
a favor sobre placas verticales, cilindros o esferas (figuras a, b y C). 


(Nu -— Nuo Y = (Nu; — Nuo)? +(Nun —Nuo)? 


donde Nu/=0, 0,3 y 2 para placas, cilindros y esferas, respectivamente (estos 


valores de Nuose conservan en el resto de las correlaciones donde 
aparezcan). La validez de esta ecuación esta bien contrastada para flujos 
laminares pero no tanto para flujos turbulentos. 


2. Flujo de capa límite con convección natural y forzada en oposición (figuras d, e 
y 1): 


(Nu— Nuo Y? = (Nu; — Nuo Y? —- (Nu, — Nuo)? 


3. Flujo forzado de capa límite sobre una placa horizontal con convección natural 
y forzada compuestas transversalmente (figuras g y h): 
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7/2 7/2 
Nu =Nuf + Nun ; L¿<L> 


? 
ES A 


7/2 7/2 12 7/2 7/2 7/2 
Nu. =Nuf +Nu» Nu  =Nuf —Nun 


donde el signo positivo se aplica cuando se calienta el fluido del lado superior o 
se enfría del lado inferior, y viceversa para el signo negativo. El límite de 


validez de esta ecuación es Nu>0. 


4. Flujo cruzado alrededor de un cilindro o una esfera con convección natural y 
forzada compuestas transversalmente (figuras i, j y k): 


(Nu-— Nuo)* = (Nu; —Nuo)* +(Nun — Nuo)! 


La raya superior sobre los Nusselt (Nu) en las correlaciones, indica que solo son 
aplicables sobre los valores medios de Nu y nunca sobre valores locales. 
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Es claro que estas reglas de combinación para la convección mixta son más bien 
burdas, sin embargo, resultan adecuadas para la mayoría de los propósitos del 
ingeniero. Las correlaciones para la convección mixta en flujos internos son más 
complicadas, para lo cual hay que acudir a bibliografía específica. Nosotros, 
simplemente no las consideraremos en este curso de introducción a la transferencia 
de calor. 


El criterio para establecer si en un proceso convectivo mixto (existen efectos naturales 
y externos) domina la convección forzada, la convección natural o es verdaderamente 
mixto, se recoge en la tabla siguiente: 


Pr<1 (metales líquidos y aire) Pr>w (aceites) 
Convección forzada Gr/Re? <<1 Gr/Pr'*Re? <<1 
Convección mixta Gr/Re? =1 Gr/Pr'* Re? =1 
Convección natural Gr/Re? >>1 Gr/Pr* Re? >>1 


donde Gr es el número de Grashof que se obtiene a partir del número de Rayleigh 
(Ra) y de Prandtl (Pr): 


Ra _ f:AT-g-L' _BATgLE, 


a r 
Pr v va 


La importancia de la tabla anterior es doble: 

- Permite establecer si en un proceso de transferencia de calor por 
convección mixta, domina la convección natural, la forzada o hay 
contribución significativa de ambas. 

- — Permite el cálculo simplificado de la convección como forzada o natural 
cuando la otra es despreciable. 


Se considerará: 


- — Gr/Re*<<1 si Gr/Re? <0,1 
-  Gr/Re? >>1 si Gr/Re? >10 


Idem para Gr/Pr'% Re? . 
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